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Abstract 

Linear acoustic wave-splitting is an often used tool in describing sound-wave 
propagation through earth's subsurface. Earth's subsurface is in general 
anisotropic due to the presence of water-filled porous rocks. Due to the 
complexity and the implicitness of the wave-splitting solutions in anisotropic 
media, wave-splitting in seismic experiments is often modeled as isotropic. 
With the present paper, we have derived a simple wave-splitting procedure for 
an instantaneously reacting anisotropic media that includes spatial variation 
in depth, yielding both a traditional (approximate) and a 'true amplitude' 
wave-field decomposition. One of the main advantages of the method pre- 
sented here is that it gives an explicit asymptotic representation of the linear 
acoustic-admittance operator to all orders of smoothness for the smooth, 
positive definite anisotropic material parameters considered here. Once the 
admittance operator is known we obtain an explicit asymptotic wave-splitting 
solution. 



1 Introduction 

The present paper derives an explicit asymptotic representation of the linear 
acoustic-admittance operator in an instantaneously reacting anisotropic me- 
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dia. This solution enables us to obtain an explicit asymptotic representation 
of the wave-splitting operators in such anisotropic media. 

Wave-splitting, or wave-field decomposition, is a tool to decompose the 
wave-field into 'up'- and 'down'-going wave field constituents in configura- 
tions with a certain directionality [19], [11] [18] [21], as e.g., a seismic exper- 
iment for probing earth's subsurface e.g., [23]. The wave-splitting proce- 
dure results in two one-way equations for the wave-field constituents. Wave- 
splitting has been used to model and analyze wave propagation in both in- 
verse problems and migration models. The method of wave-splitting has a 
long history with a wide area of applications; an overview of some of the 
history is given in [SJ. For the isotropic case wave-splitting has been used 
extensively to construct fast propagation methods [TJ] [7J [271 ES] • 

Recently, there has been an interest in methods that are almost fre- 
quency independent in calculation complexity based on wave-splitting [291 
125] . Wave-splitting has also given raise to algorithms to reconstruct material 
parameters, for example the generalized Bremmer coupling series approach 
and the downward continuation approach [TT] [23] [3U] see also [33]. Another 
application area of wave-splitting is in the context of boundary conditions 
and time- reversal mirrors [TTJ [6] . Wave-splitting methods have been imple- 
mented in several physically different contexts and for a range of constitutive 
relations: wave-splitting for wave equations [321 [33]. The electromagnetic 
equations are wave-field decomposed both for isotropic [22] E] [23] , anisotropic 
lossless (the spectral theoretical approach) [TH] and wave-splitting has been 
extended to the homogeneous lossless stratified bi-anisotropic case [26] [20] . 
Wave-splitting methods have been applied to linear-elastodynamic equations 
as for propagation on beams see e.g., [15] as well as in the half-space in 
homogeneous stratified anisotropic media [9] and up-/down symmetric me- 
dia [13]. 

One limitation to the present methods of wave-splitting is that it has 
been almost exclusively limited to isotropic media see e.g., JT9] [31] [2] [8] [23]. 
The underlying reason for this is that the wave-splitting procedures can be 
reduced to solving a certain key equation for the linear acoustic-admittance. 
This equation is almost trivial to solve in the case of isotropic media where 
it reduces to a square root of a certain elliptic operator. However, in the 
heterogeneous anisotropic media, the equation for the acoustic admittance 
is a non-linear equation, an operator Riccati equation, and this equation 
has largely resisted explicit solutions. It has been shown in [13] that wave- 
splitting methods for the exceptional case of up-/down symmetric anisotropic 
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materials resemble the methods for the isotropic case. 

The method developed in the present paper solves the above mentioned 
equation for the acoustic admittance operator in an inherently heterogeneous 
anisotropic and instantaneously reacting media, at the same time the method 
yields explicit asymptotic solutions. We expect that with the here developed 
methods, wave-splitting in anisotropic media can provide a starting point for 
wave-splitting applications in anisotropic media. We also expect that some 
of the isotropic media applications can be extended to anisotropic media. 

Before this paper there was essentially one other method for wave-splitting 
that has been extended or developed for the anisotropic case: The spectral 
theoretical approach [HJ US]- This method is a constructive method, but it 
rests on a certain spectral projection of the so called linear acoustic systems 
matrix, which can be hard both to evaluate numerically and to extend the 
method to more general constitutive relations. The leading order term of the 
linear acoustic admittance operator was obtained by [18], for the case of a 
symmetric heterogeneous anisotropic instantaneously reacting media. Higher 
order terms did not appear in their work, and seems to be hard to obtain by 
the method presented there. 

An often occurring approximation when working with wave-splitting is to 
ignore a certain lower order 'vertical' variation in the so called decomposition 
operator, see e.g., [IT] . It is known that this lower order variation can be 



accounted for in an exact decomposition for the isotropic media [321 EH] , the 
latter introduced the notation of 'true amplitude one-way wave equations'. 
An alternative approach is to include the correction term in the sources of 
the problem [111 El- m the present paper we show that the explicit asymp- 
totic admittance operator can include or ignore the vertical variation of the 
decomposition operator with minimal changes to the solution. We give both 
the solutions for both these cases in heterogeneous anisotropic media. 

This paper is organized in seven sections and an appendix, where Sec- 
tion [2] contains the explicit limitations to the material parameters of the 
linear acoustic equation considered here, as well as a reformulation of the 
linear acoustics into a form more suitable for wave-splitting. The wave- 
splitting problem is posed in Section [31 and reformulated into find a solution 
of an operator Ricatti equation for the acoustic admittance. The section ends 
with a brief comparison with earlier wave-splitting methods and their limita- 
tions. Before proceeding with the splitting-procedure, we recall the basics of 
pseudodifferential operators with parameters in Section HI The key element 
in this section is, apart from the introduction of symbols, the asymptotic 
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composition formula of symbols. 

The explicit asymptotic pseudodifferential solution to the operator Ric- 
cati equation is derived in Section [51 The asymptotic acoustic-admittance 
solution is given to all orders of smoothness in terms of a recursive formula. 
The lowest order terms are given explicitly. The result is compared with 
the known result for special cases, and our result is in accordance with and 
generalizes these known results. The asymptotic expression for the solution 
expressed in the symbol of the acoustic admittance operator is presented 
in Section [6] together with certain smoothness claims. Section [7] contains 
conclusions and reflections on our result. In an appendix we have included 
a discussion on the underlying function spaces used in the wave-splitting 
procedure as well and the normalization freedom of wave-splitting solutions. 



2 Linear acoustics and the systems matrix 

The motion of sound waves propagation through earth's subsurface is ap- 
proximated by the linear acoustic equations [211 [TJ HJ [10] . The linear acoustic 
equations are a linear system of equations in time and space which describe 
the spatial and temporal changes of the pressure, p, and the particle velocity, 
v = (v\,V2,v 3 ). In the present paper we consider the linear acoustic equa- 
tions under a time-Laplace transform, yielding the equations of motion in 
the form: 

sk(x)p(x, s) + djVj(x, s) = q(x, s), (1) 
sp jk (x)v k (x, s) + djp(x, s) = fj(x, s), j = 1, 2, 3, (2) 

where we have used the summation notation over repeated index j, k e 
{1, 2, 3}. This kind of summation notation over repeated j, k indices is used 
throughout the paper. Here x = (xi,x 2 ,Xs) is a point in space, d k ■= 
and s G C is the one-sided Laplace transform coordinate dual to time. Zero 
initial conditions have been assumed as to make — > s. Causality of the 
motion is taken into account by requiring that the Laplace-domain quantities 
are bounded functions of position in all of space for all s such that Res > 0. 
In addition we assume that 

|args|<-. (3) 

Furthermore, f k is the volume source density of force, and q is the volume 
source density of injection rate. 
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The scalar compressibility, k, and the volume density mass tensor, p, are 
the material coefficients in the equations. Sound waves in earth's subsur- 
face propagate through, e.g., water-filled porous rocks, hence p is assumed 
to be a filled 3x3 tensor or equivalently that the medium is assumed to be 
anisotropic. Both material parameters are in the linear acoustic approxima- 
tion assumed to be instantaneously reacting i.e., independent of s, and we 
assume that they satisfy the inequalities 

< k <k(x) < «i < oo, V x G M 3 

o < poICI 2 </;,/,(- < pi ICI 2 <oo,V(ei 3 ,xGi 3 . (l] 

Hence, p and its inverse are positive definite, a := p~ l , and a 33 (x) > 0, which 
is used repeatedly throughout the paper. In the upcoming calculations we 
use a rather then p in (jSJ), that is we reformulate (T5]) as 

svj(x, s) + a jk d k p(x, s) = a jk f k (x, s), j = 1, 2, 3. (5) 

The last assumption on k and a (or p) is that they are in C°° and that all 
derivatives of k, a are bounded functions. This last assumption simplify the 
upcoming microlocal analysis. 

We single out depth as the preferred direction of propagation, as for e.g., a 
seismic experiment on probing a subsurface, and choose a coordinate system 
where the X3-axis is parallel with the depth direction, also called the 'vertical'- 
axis. The wave-splitting procedure decomposes the wave-field into 'up'- and 
'down'-going wave- field constituents with respect to the vertical axis. We 
reformulate ([1]) and (jSJ) towards a form where we can apply this up/down 
splitting of the field. This rewriting has been given in e.g., [Til 0, [T5] . 
but since it is a short derivation we include it here for completeness. The 
equations are separated into two parts; the first part consists of (jSJ) with 
j = 1, 2 describing how the 'horizontal' or 'transverse' velocity components 
V\,V2 depend on pressure. The second part consists of the remaining two 
equations in (jTJ) and 

nsp + <9 3 V3 - s^d^a^dkp) =q - s _1 9 M (a Atfc / & ), 
sv 3 + a 3k d Xk p = a 3k f k , 

where we have used fl5]) with j e {1, 2} to remove all occurrences of V\, v-i in 
(HD and that Res > 0. 

We use the repeated indices p, v to indicate summation over 1, 2, just as 
we use the repeated indices j ,k to indicate summation over 1, 2, 3 throughout 
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the paper. Note in particular that the term d^a^dkP) contain derivatives 
with respect to x%. An equivalent matrix formulation of (jH]) is 



M + A)F={* S ~> ), (7) 
with F = (v3,p) T and 

\0 "33 / ' V s Ois^df, J ' 

The desired form of (j7j) suitable for wave-splitting is obtained by freeing 
c?3 from T . This is achieved by observing that the matrix operator T is 
invertible since ^33 7^ 0. We find 

(d 3 + A)F = N, (9) 

with iV = T _1 (g — s^ 1 d fl (a^kfk), «3fc/fc) T , and the acoustic systems matrix, 
.A, is given by 

— 1 i _ / ^("^"aaO SK ~ s~ x dv(Qiivdv')\ _ ( Axi A X2 
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a = 1 a = ^^^_7 3 7 _ f^"" m = Y , (10) 

in which Q^ u = a^ — a^a^a^, p,,u G {1,2}. The 2x2 matrix Q is positive 
definite since the 3x3 matrix a is positive definite. Infact, Q is the Schur 
complement of 033 in a, see e.g., [14J. An explicit proof for a similar case is 
given below in Equations (1371) and (1381) . See also [18J. 



3 Acoustic wave-splitting 

In this section we state the wave-splitting problem for the equation (Q, and 
discuss some existing solution methods and their limitations. We further- 
more reformulate the wave-splitting problem into the problem of solving an 
operator Riccati equation. 

Consider a linear invertible composition operator C = (£ + ,£~), a 2 x 2 
matrix of operators, where is the 2x1 columns of operators, which maps 
the up- and down-ward components, W = (u + ,U-) T , to the wave-field, by0 

F = CW. (11) 

1 Underlying spaces and the normalization freedom of the solution is discussed in the 
appendix. 
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We reformulate in terms of W to find 

C(d 3 +Q)W = N-(1- v)(d3C)W, (12) 
where we require £ and Q to satisfy 

v (d 3 c) + Ac = cg, g= ^ JM (13) 

for some scalar operators Q ± . Here t] takes the value of either zero or one. 
The term d 3 C in ffl2l and (Tl~3]) is below shown to be a lower order correction 
term as compared to the leading order of LQ, e.g., the operator is smoothing 
as compared with LQ. Since d$C is smoothing, it is common to ignore it, e.g., 
rj = 0, in the wave-splitting procedure. It is accounted for in this case through 
a correction term in the sources. An alternative approach is to consider the 
case with rj = 1 and adjust C with respect to these perturbations. Below we 
find solutions C v , Q v of ({TBI for both cases. 

There is a family of solutions C v which solves ( TTBl . these solutions are 
related to each other by a 'normalization' of C^. A discussion of this normal- 
ization freedom can be found in the appendix. 

To find a solution to (1131) we make the ansatz that the columns of the 
composition operator can be written as = (3^, 1) T where y^ is the 
acoustic admittance operators, it is also called the acoustic generalization 
of the Dirichlet-to-Neumann operator. Inserting this ansatz into ffl3|) and 
eliminating the occurrences of Q^r we obtain the following operator Riccati 
equation 

y±A 21 y± + y±A 22 - A n y± - a 12 - r)d 3 y± = o, v e {o, 1} (14) 

together with the relation for as 

^=^21^+^22- (15) 

Solving the equation (TTJJ yields a pair of solutions which are used to find the 
solution to the wave-splitting problem, i.e., once y^ is determined we find 
C v , Q v such that (|T3|) is satisfied. 

Before proceeding to find an asymptotic solution to fTH|) let us shortly 
discuss the existing methods for solving the wave-splitting problem, i.e., to 
find C v , Q n such that (fT3"j) is satisfied for special cases of A. For r] = and 
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the cases when An and A22 vanish is the most well known case and we find 
solutions to <HM as 

^t 0l iso = ±A 21 \A 21 A 12 ) 1/2 = ±s- 1 a 33 (a^(s 2 K - d,Q, u d v )) 1 ' 2 , (16) 
resulting in 

GU*o = ^21^t ,iso = ±K 3 V* - d,Q,M) 1 ' 2 , (17) 

since A22 = 0. This case includes the isotropic media case, where aijk{x) = 
Piso ( x )$jk an d the up/down symmetric case [13] i.e., a 3fl = = for [i =1, 
2. Note that the square root of A21A12 is a square root of an elliptic operator 
with parameter, since Q is a positive definite matrix, k > 0, | args| < 7r/2, 
and Re s > 0. 

For the case of homogeneous material parameters, e.g., k and a inde- 
pendent of x one can find solutions to the wave-splitting problem through 
diagonalization of the A in spatial Fourier domain, in this case d 3 C v = 0. 

A spectral theoretical approach to wave-splitting for r] = resulting in 
a splitting-matrix, B, has been considered by [5] for stratified layers of ho- 
mogeneous anisotropic media with instantaneous reaction, it was extended 
into heterogeneous anisotropic instantaneously reacting media in [IB] . This 
rather complex method yields, in a constructive way, the existence of the 
wave-splitting. The method constructs the wave-splitting through essentially 
a spectral projector of the acoustic systems matrix. In the linear acous- 
tics, [12] found the leading order term of the splitting operator, B, and hence 
of the admittance operator. A solution, is obtained through the equation 

yt = B 2 - 1 1 (±i-B 22 ) (18) 

where B = ^ f iR {A — A) -1 dA is a 2 x 2 matrix of operators with elements 
denoted by n, v =1, 2. The existence of the splitting matrix was shown 
under the additional requirement a[x) is a symmetric positive definite matrix. 
There are two main restrictions of the above outlined method. The first and 
more serious restriction is that the splitting matrix is complex to calculate 
explicitly, in the sense that it is derived from a spectral projection of a non- 
self adjoint operator. The second restriction is the requirement that a is 
symmetric limiting the range of material parameters that can be considered. 

Splitting for rj = 1 has to some extent been considered in frequency 
domain see e.g., [36] but it is most extensively studied in the so called time- 
domain wave-splitting procedures see e.g., [8] and references therein. Both 
approaches consider isotropic media. 
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4 Pseudodifferential preliminaries 



Our solution to the wave-splitting problem is based on the theory of microlo- 
cal analysis, in this case the theory of pseudodifferential operators. Before 
we enter into the details of finding such a solution, we need some preliminary 
notation and basic facts about pseudodifferential operators with parameters. 
We follow the notations of [28J. 

A point x = (xi, x 2 , X3) in space will be represented by a transverse part 
x± and a vertical part x 3 , e.g., x = (x±,Xs). Let the transverse spatial 
Fourier transform be given as 

u(£,x 3 ,s)= e'^' x± u(x ± ,x 3 ,s)dx ± , (19) 

where £ = (£1,^2) is the Fourier dual to x±. Since A^ u for all \i ,v =1, 2 are 
differential operators with smooth coefficients, we find the left-symbol, a 12 , 
of A12 by a 12 (x,£; s) = e~^' x± Ai 2 (x, di, d 2 , s)e^' x± . Resulting in 

a 12 = sk + s^Q^^ - s~ 1 {d lx Q lxv )ii v . (20) 

Similarly 

On = i^a^a^ + (d^a^a^)), a 21 = sa^, a 22 = i^a 3fl a^. (21) 

The above symbols are symbols with a parameter, s 6 C. Since s corresponds 
to the time derivative, it is necessary to consider the parameter s as if it is 
of the same order as £. 

The symbols are decomposed into order of homogeneity in s, £. To illus- 
trate this concept, we separate a\ 2 into its poly-homogeneous components in 
s, £: a\2 = 0,12-1 + «i2;0! an d similarly for the other components of a. The 
explicit form for each of the poly-homogeneous terms of a is 

011,1 = i^aA"3 a 33 a n,o = ^(a^agg 1 ) (22) 

022,1 = i£,n a 3/J> a 33 °21,1 = Sa 33 ( 23 ) 

012,1 = sk + s~ l Q^^ u , 012,0 = -s^d^Q^i^. (24) 

The composition formula for symbols see e.g., [28] is formulated as follows: 
Let P and Q be two (properly supported) pseudodifferential operators in M? 
with corresponding symbols q{x±,£), p(x±,£) respectively. The composition 
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formula yielding the symbol, r, corresponding to R = PQ, is given by the 
asymptotic relation 1.3.6], [3J 1.1.5] 



(25) 



where (3 is the multi-index (3 6 N 2 . Here /3 = (Pi,f3 2 ), ft = and 

The correspondence between a symbol, r, and its operator, i?, acting 
upon some function it is 



(2*) s 



e^<r(x ± ,Om de 



(26) 



Note that partial differential operators with smooth coefficients are properly 
supported, and by a uniformity argument this property extend to partial 
differential operators with parameters, it is also true for each given branch 
of the square root of elliptical operators with parameters. 



5 Asymptotics of the acoustic admittance op- 
erator 

In this section we derive the asymptotic expansion of the admittance operator 
y. This is rather long recursive expressions, where the next term is expressed 
in terms of all earlier terms. Once we derived the general expression we 
compare the derived result with existing results for a pair of particular cases. 
Henceforth we suppress the index r\ for readability. 
Starting from the operator Ricatti equation 

y(A 2 iy + A22) - A n y - a 12 - v (d 3 y) = 0,^0,1 (27) 

we are interested in deriving the asymptotic expansion of the symbol y corre- 
sponding to a solution y of (12"7|) . We expect to get two (separate) solutions 
y + and y~ from (1271) . The normalization freedom of £ imply a normaliza- 
tion freedom of y and hence we expect an infinite number of solutions in two 
equivalence classes, see Appendix |A] However, for our purpose it is sufficient 
to find one pair of solutions to fl27|) . e.g., y ± . 

To translate (}2"T]) into an equation for symbols, we repeatedly use the 
composition relation (|25[) . Let's denote the symbol corresponding to y by y. 
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We start with the middle part of (12"7|) : IZi := —Any — A12, its symbol, r 1; 
has the asymptotic expansion 

n ~ -i^oc^a^y - d^a^a^y) -sn- s^Q^^ - s^d^Q^i^ (28) 

after simplification. The simple form is due to that An is a first order 
differential operator. Similarly we study the first part of ff27|) and denote this 
term IZ2 '■— y (A21 y + A22), its corresponding symbol, r 2 , can be obtained 
directly from fl25|) . Once again due to the simple form of A21, we find that 

Solving the operator Riccati equation on a symbol level is equivalent to 
solving the equation 

r\ + r 2 - 77ft.?/ = 0, (30) 

or equivalently solving 

-w^ d ty)(- d xfi.sa^ y + i^a 3fl a^) - i^a^a^y - d^a^a^y) 
p p - 1 

- - s _1 <5^^^ + s^d^Q^y - r]d 3 y ~ (31) 

with respect to y. We observe that the ansatz y = y + y„i + y_ 2 + . . ., where 
each term, y_ n is poly-homogeneous in £, s of order — n, yields a consistent 
solution of fl3Tl) . With this observation, it is also clear that the term d$y is 
a lower order term, see also Section [61 

Inserting the expansion of y into (I3"TT) enable us to solve the equation in 
an iterative manner, first we collect all l:st order terms: 

saisJ/o + ^"asC^ ~ "^3)2/0 ~ sk + s^Q^i^ (32) 
with solutions 



Vo ~ s 1 ( - - "^) ± 7i) • 



(33) 



where 

7i:=a^V« + <UU»0 1/a (34) 
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with 

Qnvin^v = {Q^ - ^( a 3 M - a^)a^( a ^ - Oiu3))^u (35) 

= (<V ~ \( a ^ + a 3^) a 33( a f3 + CX3u))€n€v (36) 

Note that Q^ v is positive definite, this follows from the assumption that a is 
positive definite and the observation that for any vector £ e R 2 we can find 
a vector ( e R 3 such that 

Qfivitiiv = 0!jk(jCk (37) 

and since a > we find that Q > 0. The relation between £ and £ is 

C = -^aasCofy + "ais)^)- ( 38 ) 

This proof is similar to the proof that if a > then Q > given in [TH] . 
The zero order terms in fl3TD are: 



«33 1 2/-l( 2s 2/0 = + i ^(«3 At -tt At 3)) + Yl ^[(^0 = )(y^) /3 ( Sa 33 1 2/0 = + i ^ a 3 At a33 1 ) 



- d^a^a^y^) - VdaVa s ^Q^v (39) 

Once we observe that 2sy$ + ^(a^ — a^ 3 ) = ±27! and recall that 71 is the 
first order symbol of a square root of a second order elliptic operator with 
parameters since |args| < it/2 and Res > and hence invertible [28]. We 
find that the next order term, y_i has the asymptotic behavior: 



33 / 



(40) 



We proceed similarly for each poly-homogeneous order in (131 j) and can 
thus from the — n:th order terms find the y^ n _i solution for n > 0. The 



2 Recall that we always sum over 1, 2, 3 for repeated the repeated indices j, k, and over 
1, 2 for the repeated indices /1, 
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solutions are expressed in terms of y_ n , . . . . Its explicit expression is 



k=l \/3\=k j+m=k—n—l 
—n<j,m<0 

With the results in (l33l) . (|4T)]) and (|4T!) we have thus obtained an asymptotic 
representation of y^ and hence of y 1 * 1 . This solution solves the operator 
Riccati equation (B71) . 

Comparing with the case of a symmetric a = a(x) as detailed in [TB] we 
find that their leading order term agrees with our leading order term for the 
case when r/ = 0. Furthermore, for the simpler isotropic case the above result 
agrees with the leading order term reported in [23] . 




6 Claim of lower order 

In this section we use the above derived acoustic admittance symbol to ex- 
press the respective (matrix) symbols of the wave-splitting solution. We also 
show the claim that <9 3 £ is smoother than CQ. 

The above introduce ansatz = (y ± , I) T does indeed give us a solution 
to the splitting problem. Starting from (fT5l) we find that the diagonal terms 
in Q have the symbol, g ± with poly- homogeneous expansion 



= 9i+9o+9-i + --- } (42) 



9- 1 = asiisVo + i^«3 M )> 9-n = ®33 S y-n-l, n = 0, 1, 2, . . . (43) 



Hence g ± is expressed as a poly-homogeneous sum with leading order term 
of poly-homogeneous order one. 

Similarly we find that the composition operator £ with symbol, £, of the 
form: 




j+k=— n— 1 
— n<j,k< — 1 




(41) 




where 




(44) 
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where 



I y ly *-n=[o o"J' n = > 1 > 2 '~ ( 45 ) 

We claimed in Section [3] that the term d 3 C is of lower order than CQ, we 
are now in the position to show this claim. We do this by explicitly giving 
the symbols of d 3 C and CQ. The symbol of d 3 C can be written as 

where 

d 3 y ± ~ s _1 ( - ^i^<9 3 («3^ - a M 3) ± tt-(s 2 33 k + (d 3 Q^)^u 

+ terms of order -1 and lower. (47) 

The symbol, p, of P := is 

P ~ E(^)(i^)^ ~ + l-i9i + to9o + Yl ^l ){\d x f 9l 
P 1 W=i 1 

+ terms of order -1 and lower. (48) 

where £ m is the 2x2 matrix see (jUJ) and g m = diag(g+, g~), with ^ for 
m G 1, 0, —1, —2, • • • as given in (142]) . 

Comparing the result for the respective leading orders, of d 3 £ and p, we 
find that 8 3 £q is a matrix containing the terms (d 3 y^,0) while the leading 
expression for p is £ gi with terms {y$gf,gf). We find the claim to be 
shown, CQ has leading order symbol of order one whereas d 3 C has leading 
order symbol of order zero in the poly- homogeneous expansion of the symbols 
in s, £. Hence, d 3 C is a smother term as compared with CQ as claimed. 



7 Conclusion 

The main advantage with the above procedure to derive as compared 
to the spectral theoretical approach is that it gives an explicit asymptotic 
representation of y ± through its asymptotic symbol expansion y ~ yo+2/-i + 
. . . as given above. The method yields, as far as the authors know, the first 
explicit form of y_ n for n > in anisotropic media. 
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The presented method is not limited to material parameters for which 
certain spectral properties of A exist, in our method it suffices to consider 
material parameters k(x), at(x) such that 033 and s 2 k + d x Q^d^ have well 
defined inverses and square roots. 

The flexibility of the presented approach is shown by the possibility to 
include the often neglected term due to the vertical derivatives of the com- 
position operator with minimal changes to the derivation procedure. The 
simplicity of the method indicates that it might be extended to the case of 
the electromagnetic equations, possibly with anisotropic conductivity, as well 
as to linear elasticity. 

A Non-uniqueness of wave- splitting solutions 

To formalize our discussion of the normalization freedom of wave-splitting 
solutions we introduce the function space H n := (H n (M 3 , C), H n+1 (M 3 , C)) T , 
where H n is the n:th order Sobolev space of square integrable functions, 
net, see e.g., [3]. 

The normalization freedom implicitly appear in the first step in phrasing 
the wave-splitting problem: 

F = CW. (49) 

Here we need to specify to what spaces F and W belongs. Let's start with 
the case outlined in the paper. We have implicitly assumed that the sources 
iV in ([9]) belong to H n_1 , for some n6l and hence, due to the derivatives 
in A, we find that F e H n . The obtained £ with symbol given in (|44|) maps 
H n+1 to H n , and the obtained wave- field constituents satisfy W G H n+1 . 
Indeed if we use the notation C n to indicate that C n : H n \— > H" _1 and 
similarly for Q we can rephrase the requirements on C, Q as: find £, Q such 
that 

e {0, 1} (50) 

with Q n+ i diagonal. 

To make the normalization freedom apparent we replace (j49p with F = 
C n+ iW, where 

l m = C m M~\ W = M n+ iW, m,neR (51) 

for any invertible operator M m from H m to some desired space W m for some 
m£l. 
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The natural question arises: what are the requirements on Af so that 
C and its corresponding Q is a wave- splitting solution to (ITT1) - (IT3"1) . These 
operators M represents the normalization freedom of the solution of the wave- 
splitting problem. It is rather straight forward to implicitly characterize 
which M that provides a wave-splitting solution, e.g., we substitute C,G in 
(JE0 into (JI2P and (JX3J) to find 



where we have identified Gn+i '■— AfnGn+iNn+i ~ r){dzNn+i)N n +x- The re- 
quirement on Af m is that Gn+i remains diagonal for diagonal Gn+i- In ad- 
dition as an implicit requirement on A/", for W G W n+ i we require that 
(dzNn+ijN'^W G W n , as indicated in the restriction of the C n+ i term to 

l n in front of (9 3 A^ + i)AC+i in §2]). 

The normalization freedom above is rather large, there are several diag- 
onal (and anti-diagonal) M which satisfy the above conditions, e.g., M = 
diag(m, m'), m, m' G M and M = diag((y + ) p , iy~) p ), p, p' G R. The nor- 
malization freedom was used in earlier wave-splitting papers, for the isotropic 
case [11] utilize the normalization freedom to derive wave-splitting solutions 
with self-adjoint admittance operator. Normalization was also discussed in 
the anisotropic case [18] where the spaces H n also appeared. 

One case of some interest is the normalization that switch from the 
acoustic-admittance representation to the acoustic-impedance representa- 
tion. Here we are interested in a normalization operator 



and W G H™. Thus Af n : H" i— > H n 1 . The explicit form of the symbol 
of y ± ensure the implicit assumption that {dzN' n+ \)N~+- i W G H n and that 
C n+1 : H n ^ H n . We obtain 



(r](d 3 C n+1 )+AC n+1 )W 



Cn [KGn+lK+l - VidsMn+lWn+l] W = CjW , (52) 



(53) 




(54) 



and 



G 
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